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A NOTE ON 

SOUND RADIATION INTO A UNIFORMLY FLOWING FLUID 

by 

Harold Levine 


Introduction and Summary 

There exists a considerable body of analysis pertaining to the 
sound radiation generated by mechanical or vibratory source arrangements 
and, in some instances of planar or piston types, the results encompass 
all magnitudes of the wavelength or frequency as well as disposition of 
the observation site. Fewer details are available in the circumstance 
of relative motion between the source and its surroundings and, partic- 
ularly, with an aptness for both the compact and noncompact characteriza- 
tions of a given source model. The comparative increase in average total 
radiated power from infinitesimally small or point sources of peri- 
odic strength, on passing from rest to steady rectilinear motion, de- 
pends only on the Mach number in the latter state; and a formally analo- 
gous rise in output is linked, by Ffowcs Williams and Lovely (1975) , 
to the presence of a steady parallel flow past a rigid plane wall in 
which a compact circular piston executes normal oscillations. It is 
the intention here to widen the perspective of effects connected with 
such a background flow by regard for an elongated or strip piston, 
which prototype permits a straightforward and general analysis. 

The time average net power output of a strip piston that vibrates 
normally to itself with uniform amplitude, has a finite width 5 and is 
located in an otherwise rigid plane wall, admits different expansions 
according as 
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k6 « 1 


M < 1 


(i) 


or 


Mk6 « 1, 


( 11 ) 


or 


kfi » 1, M < 1 


(iii) 


respectively, where k,M are the wave and Mach number, the last being a 
measure of the ratio between flow and sound speeds. In the absence of 
flow, M vanishes and classical (or linear acoustical) theory predicts 
that the total rate of energy radiation (which is distributed nonuniformly 
over a semicircular range) amounts to 

P = p^kc^(k5) j / J^(v)dv - J^(k6)j (1) 

1 ? ( " J- (V) I 

= ip kc^(k5) |l - / - i - - dV (2) 

^ ° k6 ^ ) 


per unit amplitude and length of the piston; here p^ specifies the equi- 
librium density of the medium, c is the sound speed and J^, are 

Bessel functions with the designated orders. The development which 
suits a compact piston. 


P = i p kc^(k6)^ 
4 ^o 


(k6)- 

24 


+ 0(k6) 


k6 0 


(3) 


IS a ready consequence of (1) , while that obtained from (2) , 

P = y pQkc^(k6) jl - (k(5) cos^k6 - J) + 0(k6) k6->™ (4) 
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represents the noncompact alternative. 

It appears, in what follows, that the development containing 
powers of the Mach number and conforming with (ii) , 


k6 


P = 4 P kc^(k5) I / J (v)dv - J. (k6) 
2 o I 1 o 1 


MkS « 1 


+ M 


1 3 ^ 

j k6 J^(k6) + I Jj^(k6) + ^ 


o(m^) , 


(5) 


extends the prior result (1) and implies that 


P = i p^kc^k6)^ 


1 + 3M^ 



( 6 ) 


in the limit k6 0; if k5 » 1, on the other hand, the expression 
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k6 00 (7) 


takes the place of (4) . The principal effect of flow in (6) is merely 

2 

to raise the long wave power estimate by a numerical factor 1 + 3M , 

whereas the comparison of (4) and (7) discloses that the flow brings 

about a significant change in the relative order of magnitude of the 

-3/2 -1 

second term in the short wave estimate, namely from (k6) ' to (k6) 

Inasmuch as the second term of (7) is, in the absolute sense, indepen- 
dent of the piston geometry (and, specifically, of its breadth) , a cor- 
responding feature may be presumed for other shapes. A further 
estimate. 


P = V P„kc^(k5)^ 


1 + |m2 


(i- 


2 

M / 


\5/2 


+ 0(k6) 


k6 « 1 


0 < M < 1 , (8) 
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which accompanies the stipulation (i ) , has a more extensive validity in 
respect to the Mach number than does (6) . 

The Analysis 

Let the undisplaced piston occupy a strip |x| < a, - < y < oo, in 

the fixed plane (or wall) z = 0, and assume the existence of a steady 
flow with uniform speed U = Me in the parallel (or x-) direction. 

Given the normal piston displacement (exclusive of a periodic time fac- 

. -10)t. 

tor e ) , 


z = n(x) = / A(C)e^^^d? 


(9) 


where 


1 r . -li^x , sin i;a 

- 2? / = —FT” 


(10) 


for the 'top-hat' profile 


n(x) = 


If |x| < a 
0 , |x| > a 


( 11 ) 


then the corresponding acoustical excitation, derivable from a velocity 
potential (})(x,z)e is implicit in the convected wave equation 


V^({) = f-ik + M (t> 

dx/ 


k = w/c 


or 


( 12 ) 


I 2] 3^(j) 9^c() 3(|) 2^ 

(l - M / — ^ ^ + 2iMk + k 4) = 0 


3x 


3z' 


3x 


z > 0 
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together with the boundary or wall condition 


3z 


z=0 


9n 

-x(on + u ^ 


( 13 ) 


and an outgoing wave condition at infinity. 

A solution of the linear partial differential equation (12) has 
the form 

00 

(f)(x,z) = / B(^) exp{iCx + iK(C)z}d? (14) 

— 00 


where 


K^(?) = (k - , arg K(?) > 0 


(15) 


and compliance with the boundary condition (13) necessitates a propor- 
tionality 


B<0 


(16) 


between the weighting factors A (^) , B(C) of the respective integrals (9), 
(14) . The normal derivative of 4>, as computed from (14) , (15) and (10) , 
proves to be 


(JU 

III = i / B(?X(?)e^^^dC = - i (w - CU) A(C) e^^^ d? 


Z=0 - <X> 


r sin^ i^x lu J 

IT r IT J 


^ {sgn(a + x) + sgn(a - x) } + u{6 (x + a) - 6(x - a) } (17) 
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where sgn V = ± 1, V ^ 0 and 6 signifies a Dirac delta function; it is 
thus evident that 3(|)/9z vanishes on the rigid section of wall (|x| > a) , 
has a uniform magnitude over the piston (|x| < a) and manifests the 
joint presence of a line sink/source at the leading/trailing edges of 
the piston (|x| =:fa). 

The time average power output of the piston is conveniently found 
through its rate of local working on the adjacent medium (with an equi- 
librium density p^) , viz. 

P = ^ p w Im / 4>*(x,o) dx (18) 

2 ioo 9z 


and, consequent to the use of the representation (14) for 4» and the con- 
nection (16) , this becomes 

“ 2 

P = irp^u Re / |b( 0 I K(C)d? 

—00 


2 f 

= up to(kc) Re I 

O — GO 




|(k - 


1/2 


A(i;) 


dx 


(19) 


the same result obtains if the transformation of (18) commences with a 
version 

OO 00 

p = -i-pwim fdx / B*(?)e"^^^ I “ [sgn (a + x) + sgn (a - x)] 

2 o -00 ( 2 ■■ 


+ u[6(x + a) 



dC 
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that incorporates the explicit determination of 3(j)/8z given in (17) . On 
substituting the expression (10) for A(C) and simplifying (19) , it ap- 
pears that 


P = 


Poke 


3 i-M 


(1 - MX)' 


Sin Oeax) 


^ , 1/2 
--i- {(1 - MX) ^ - X^} 

1-M 


dx 


Pokc- 


(1 - Mg)^ 


sin^/lca 

\ 1 - M^/ 


7T 


/l - ^ 


/TT 




M < 1 


(? - M) 


2 

in which the change of variable X = (5 - M)/(l - M ) figures. 

To reduce the latter integral put ^ = cos 9 and introduce the ab- 
breviations 


a = ka , 



( 20 ) 


whence 


P = ^ 




(a,M) - 2MFj^(a,M) + M F 2 (a 


,M)| 


( 21 ) 


with 


„ , 'f n^ sin^ (otB(cos 0 - M) ) 

F^(a,M) = J cos 0 

o (cos 0 - M) 


d0, n = 0,1,2 


( 22 ) 


An alternative single integral representation for each of the functions 

F (a,M) , which IS the more suitable in seeking estimates appropriate to 
n 

different hypotheses regarding the magnitude of the argument variables, 

rests on the fact that the second order a - derivatives of the F are 

n 

explicitly known in terms of trigonometric and Bessel functions, viz. 
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(23) 


and 


3^F 


3a 


2tt 3^ cos(2aeM) J^(2a3) 


3^f 


3a 


i= 2 tt 3^ sin(2aBM) J^(2a3) 


(24) 


3^F, 


3a" 


= 2f 3^ cos(2a3M) |j^(2a3) - Ji(2a3)| 


(25) 


The two-stage integration of (23) , which serves as a model for that of 

(23) , (24) and is rendered definite by the universal conditions 
3 f^(0,M) 

F (0,M) = = 0, begins with the expression 

n 3a 

3f a 

-^ = 27t3 / cos(23mt) j^(23x) dx 

o 

and carries on via the sequential formulas 
a V 

F^(a,M) = 27 t3^ / / cos(23mx)j (23x) dx 

o ^ o 

o o 


.( “ 

a J 

= 2-rr3 ]a / cos(23mx)j ( 23 x)dx - / x cos(23mx) j ( 23 x)dx [ 

' o 

o ) 

2a3 

2a3 

= Tra3 / J (v)cosMvdv- ^ 

J VJ (V) cos MV dV 

o ° 2 

0 

2a3 1 

i ^ 2 

= ^3 / Jq (V) cos MV dV - — 

o 1 

#- J J (V)sinMVdV-4a3 Msin(2a3M)J (2a3) 
dM *^0 o 

1 o 


= TTa3P(a,M) + 2Tra3^M sin(2a3M) J (2a3) (26) 

2 oM o 
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where 


2a6 j 

P(a,M) = / j (V)cos MV dV = 3 - / J (V)cos MV dv 

O ° 2aB ° 


(27) 


and M < 1 

2a3 ~ 

Q(a,M) = / J (V) sin MV dv = - f J (V) sin MV dv . (28) 

o n o 

o 2a3 

The analogues of (26) for and F 2 are 

F^(a,M) =7Ta3MP(a,M) - jQ(a,M) -jM^Q(a,M) + 2ira3^M^ sin(2a3M) J^(2a3) 

(29) 

and 

F„(a,M) = 7Ta3M^p(a,M) - ^MQ(a,M)--^ ^Q(a,M) + + m^) sin(2a6M) j (2a3) 

Z Z Z qJ4 ^ o 

+ 7ra3 cos(2a3M) J^(2a3) + j|l - cos(2a3M) J^(2a3)l (30) 


respectively. 

Hence, the three integrals F^, F^, F^ are charapterized , apart 
from explicit terms, through the simpler pair, namely P(a,M) and Q(a,M); 
and the representation for the power which involves the latter turns out 
to be 

p = p(a,M) +3^ m(i - ^)Q(a,M) - j ^ Q(a.M) 

_3_ 3 

+ a3^ m( 2 - M^) sin(2a3M) J (2a3) + a3^ v? cos (2a3M) J, (2a3) 

O 1 

+ y 3*^ M^(l - cos(2a6M)J^(2aB)) j . (31) 
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Multiple differentiation of P(a,M) and Q(a,M) with respect to M, at the 

specific value M = 0, permits the formation of developments containing 

powers of the Mach number , which commence in the fashion 

2a i 2a \ . 

P(a,M) = / J^(v)dv + y M^|-4a^Jj^(2a) + 2aJ (2a) + / jQ(v)dv[ + o(m'^) 
o ( o ) 

and 

Q(a,M) = 2aMJ^(2a) + j M^|-4a^Jj^ (2a) + 8a^J^(2a) +4aJ^(2a)| + o(m^) . 

Employing these estimates and approximating the other terms of (31) in 
accordance with the hypothesis that a3M << 1, or simply Mka « 1, the 
development (5) cited in the introduction emerges after the dimension- 
less magnitude 


2a = 2ka = k6 


is redefined in terms of the strip breadth 5- 

The same development for the average rate of energy radiation can 
be secured, as a matter of technical interest, by the 3 udicious extrac- 
tion of finite parts for divergent integrals; thus, recalling the 

Legendre polynomial generating function 

— 00 

(l - 2zt + t^) ^ ^ Pj^(z)t"^ , |t| < 1 (32) 

n=0 


and writing 


|(k - M?)^ 
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with the identifications 


t = 




z = 




/T^ 


(33) 


it IS formally inferred from (19) that 

i2 


p^O)(kc) ' 


P = 


Re / 


- f )‘ 


Sin Ca 








(34) 


k 2 

2 , 2 r sin 1 

= -p^ 0 )(kc) / ; 




2 „ 2 , 1 

k2 " 2 ^ 2 _^ 2 \ ^ 2_^2 


- 1 




d? 


(35) 


Since the restriction 1 1 1 <1 which assures convergence of the expansion 

(32) does not hold uniformly in the context of the pertinent formulas 

(33) , (34) all save a pair of terms in (35) possess nonintegrable singu- 
larities (at ^ = k) . The proper integrals are 

k 2 ^ 2 k^ ^ 

F(ka) = J ^ ^ ka I / J (v)dv - J, (2ka)/ 

° ' /k= - ° ’ 


and 


G(ka) 



• 2 - 
sin Ca 

/k^ - 


dC 



Jo(2ka)) ; 


after noting an alternative expression 


k 

I 


o 


sin Ca 





2r 

sin Qa 
2 
C 



dC 
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for the first improper integral, its finite part may be identified as 


dk 


k / 


2 

sin Sa 


dC 


o C 


/ 7 T 7 


- - -I 

dk ( 


d /F(ka)\ TTa 


k f 2k 1 


J (2ka) 


Likewise, the prescription 


r sin ^a 1 


/ 7 T 7 \ k - 


- 1 dc 


= f s^n^ga 1 1^ _ ^ k 




k^ - 


+ 3 




d^ 


dk 


2 

k2 J s in Ca 


° cV 


r/ 


d 

dk 


k f s i^ ..-- dg: 
° cAk2 - e 


dk^ 




2 

= TTa J^(2ka) - II Jj^(2ka) 


supplies a value for the other improper integral which appears in (35) . 
When the various integral determinations are brought together with (35) , 
the resultant power estimate, 

( 2ka 

/ J (V)dV - J, (2ka) 
o o 1 

+ kaJ (2ka) + 4 <2ka) 

1 O 2 1 

Mka « 1, 


1 - 


Jo(2ka) 

__ 


+ 0 M 


agrees with (5) . 

Turning attention next to a noncompact piston, typified by the in- 
equality a = ka » 1, and having regard for the ultimate stage of the 
successive equalities 
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/ J (v)cos MV dV = / /vj (V)) 

Y y \ / 


cos MV 
V 


= - J (Y)cos MY + / 

Y I V 


cos MV . „ sin MV , 

+ M n > dV 


Y 

00 J^(v) °° 

= - J (y)cosMY+ / — n — cos MV dV + M j d/-J (V)| sin 

-*■ y ^ -V \ ' 


■} 


MV 


“ Jn (V) 


1 ' ' 2 

= - J^(y)cosMy+ / — — — cos MV dV + MJq (Y) sin MY + M / J^(V)cos MV dV 


Y 


Y 


It follows that 


( Jq (V) cos MV dV = 2 

Y 1 - M^ 


^ 00 “1^ 

J (Y) cos my + MJ (Y) sin MY + / — cos MV dV 

1 O y V 


(35) 


whence the function P(a,M) defined in (30) is expressible as 

J I ” J (V) 1 

p(a,M) =3 + 3jji(2a3)cos(23M) -MJQ(2aB)sin(2a3M) - | —^ 5 — cosMVdvj . 

I 2a3 ) 


(36) 


Correspondingly , 


°° J, (V) 


/ J (V)sinMVdV= < -Jy (Y) sin MY - MJq (Y) cos MY + / — ^ 

Y ° 1 - M^ I Y 


sin MV dV 


(37) 


and 


Q(a,M) = 3 


“ J,(v) I 

J (2a3) sin(2a3M) + MJ (2a3) cos (2a3M) - f sinMVdvf 

^ ° 2a3 ^ ’ 


(38) 


13 



from which it is deduced (on employing the result d3/dM = 2B M) that 


do 2 2 2 

^ = g M J (2ag) cos (2agM) + 3 MJ (2ag) sin(2ot3M) 
dM o 1 


2 7 


+ 2a3 MJ (2a6)sin(2a3M) + 2a3 J, (2a3)cos(2a3M) - m3 J — ; — smMVdV 
° 1 2a3 ^ 


( 39 ) 


Substituting the expressions (36), (38) and (39) for P (a,M) , Q (a,M) and 
dQ/dM into the formula (31) constitutes the last step of our power analysis, 
and yields the generally valid representations 


P/p^kc^ = a + I 3*^ 


L 

- aB’ / — 


J, (V) 


2a3 


1 ^ 

^ cosMVdv + — m 3 J^ (2a3) sin (2agM) 


1 ^ 7 

~ TT B M J sin MV dV 

2 2a3 


(40) 


and 


k6 , Q 

a = — = ka, 3 


1 - M 


2a3j^(v) ^ , 

p/p^kc'^ = a3^ / — — cosMVdv + J m 3^ Jj^(2a3)sin(2a3M) 


1 I, 

o 


sin MV dV 


(41) 


where the former lends itself to estimation if a3 » 1 and the latter if 

a3 << 1. In particular, when a 0 and M does not assume values close to 

unity (or 6 close to infinity) the representation (41) implies 

5 1+12 

P/p^kc^ = (i + 1 M^) + o(a'^) = i (k6) ^ + O(kfi)'^ 


(l - m2) 
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M O 


= j (k6)^|l + 3M^ + o(m^)|’ + 0(k6)' 


as stated previously. For small values of M, furthermore, the develop- 
ments 


2a3 J]^(V) 

/ — cos MV dV = 



2a 

= / J (v)dv - J. (2a) 

n -L 



2a > 

2aJ^(2a) + 2J^(2a) - / J^{V)dv( 



and 

2ag J (V) 2a I . 

/ sin MV dV = M / J, (V)dV + 0\M ) 

o o 

= M|i - J^(2a)[ + o(m^) 

support, in conjunction with (41) , the propriety of the result (5) . 

The terms of (40) are arranged in a sequence that befits their 

relative importance in the limit a M ^ 0; namely, the first and 

second are 0(a) and 0(1) , respectively, while the third and fourth 
-J, -3. 

are 0(a ) and the last is 0(a ■*). Of particular note, as previously 

stated, IS that the flow manifests it presence, to the leading order 
in the average power output at short wavelengths, by a contribution 
which does not involve the piston scale. Estimates for the terms of 
(40) which contain Bessel functions are readily gained through the 
utilization of large argument asymptotic forms pertinent to these 
functions. 
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